Orientation and Motion of Water Molecules at Air/Water Interface 
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Analysis of SFG vibrational spectra of OH stretching bands in four experimental configurations 
shows that orientational motion of water molecule at air/water interface is libratory within a limited 
angular range. This picture is significantly different from the previous conclusion that the interfacial 
water molecule orientation varies over a broad range within the vibrational relaxation time, the only 
direct experimental evidence for ultrafast and broad orientational motion of a liquid interface by 
Wei et al. [Phys. Rev. Lett. 86, 4799, (2001)] using single SFG experimental configuration. 

PACS numbers: 33.20.Tp, 42.65.An, 68.35.Ja, 82.30.Rs 



Due to lack of direct experimental measurement tech- 
niques, motion and dynamics of molecules at liquid in- 
terface are yet to be well understood 0. In the past 
decade, air/water interface, as one of the most impor- 
tant liquid interface, has been intensively investigated 
both theoretically and experimentally ^ S S M S 
0, H, 10, 11, 12, 13, 14). Experimental method such 
as X-Ray reflection can be used to get the liquid surface 
roughness but orientational structure and dynamics 
information at molecule level can only be obtained by 
Sum Frequency Generation (SFG) or Second Harmonic 
Generation (SHG) because of their submonolayer sensi- 
tivity and interface specificity 0. It has been gener- 
ally accepted that the interfacial water molecule has one 
free OH bond protruding out of the interface. However, 
the existence of water molecules with both OH bonds 
minting to the vapor phase is still a controversial issue 
^ m 0, Q|. The dynamics of water molecule at 
air/water interface has been discussed theoretically in lit- 
eratures The only direct experimental 
study so far concluded a fast orientational motion varies 
over a broad range of 102° centered at the surface normal 
{9m = 51°) within a time scale comparable or less than 
0.5 ps, by Wei et al. from the fact of the vanishing of the 
SFG vibrational spectra of the free OH stretching mode 
of the interfacial water molecule around 3700cm~^ in 
the sps polarization (denoting s- p-, and s- polarized sum 
frequency output, visible input, and infrared input, re- 
spectively) y . This result indicated a very dynamic and 
disordered physical picture for the air/water interface. 

Femtosecond pump-probe experimental studies have 
shown that orientational relaxation dynamics of water 
molecules in bulk water is as fast as Ips, and also the 
motion of non-hydrogen bonded water molecules is faster 
than hydrogen bonded water molecules 15j. Therefore, 
it is reasonable to predict very fast orientational mo- 
tion of interfacial water molecules, which is less hydrogen 
bonded as bulk water molecules. However, theoretical 
simulation by Chandra et al. argued for slower dynam- 
ics for interfacial water molecules than that of bulk ones 



(10.5ps vs. T.lps) 0. Simulation by Benjamin et al. 
also suggested that the lifetime of hydrogen bond at inter- 
face between water and an organic liquid is significantly 
longer than that of bulk water molecules . No matter 
how, these all suggested ultrafast orientational motion 
of interfacial water molecules. Most importantly, ultra- 
fast orientational motion may have effects on vibrational 
spectral width, as well as the strength of vibrational res- 
onance in different input/output polarization combina- 
tions, provided that the motion is faster than the time 
scale of l/Fg (the rapid motion limit) and the motion 
covers a broad range centered around the surface nor- 
mal, as reported by Wei et al. 9J. 

Recent quantitative analysis of data in SFG vibrational 
spectroscopy have suggested that vapor/liquid interface 
are generally well ordered, and sometimes even with anti- 
parallel double-layered structures 0, 0, 0, 0| . It has 
been generally accepted that liquid interface with strong 
hydrogen bonding between molecules should be well or- 
dered []|. Therefore, it is surprising to comprehend the 
picture dictated by ultrafast orientational motion in such 
a broad angular range. 

Here we show from SFG-VS measurements in four sets 
of experimental configurations that Wei et aL's analysis 
13 can only be considered valid for the particular SFG 
experimental configuration. We conclude that the ori- 
entational motion of water molecules at the air/water 
interface can only be libratory within a limited angular 
range around a tilting angle around 33°, instead of being 
around the interface normal. Therefore, the air/water in- 
terface is much better ordered than previously suggested. 

The detail of the SFG-VS experiment was described 
elsewhere [i^. Fig. ^shows SFG spectra of the air/water 
interface in reflective geometry with .ssp, ppp and sps po- 
larization combinations at four sets of incident angles for 
the visible(Vis) and IR laser beams, namely, Config.l: 
Vis== 39°, IR= 55°; Config.2: Vis= 45°, IR=55°; Con- 
fig.3: Vis= 48°, IR= 57°; Config.4: Vis= 63°, IR= 55°. 
It is clear that there is strong dependence of ppp and sps 
spectra on the incident angles, especially on the visible 
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beam. The sps spectra in Config.4 shows clearly much 
larger intensity at ~ 3700cto~'^ than those of the other 
three configurations. Even though the sps intensity is 
generally much smaller than that of ssp and ppp polar- 
izations, the sps intensity is well above the experimental 
error, and have been repeated more than a dozen times. 
It is important to see that the spectra in Config.2 re- 
sembles closely the data obtained by Wei et al. with 
Vis= 45°, IR=57° 9]. 
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FIG. 1: SFG spectra of air/water interface in different po- 
larization combination and experimental configurations. All 
spectra are normalized to the same scale. The solid lines 
are globally fitted curves with Lorentzian line shape function. 
Note the different error bars for graphs in different scales. 



The ssj) spectra of all four experimental configurations 
overlap with each other well within experimental error 
when normalized to the ~ 3700cm~^ peak intensity. The 
difference of their absolute intensities can be quantita- 
tively accounted from the Fresnel coefficients with differ- 
ent incident angles. These facts indicate that SFG data 
with different incident angles are from interfacial layers 
without detectable bulk contributions [2H . 

Quantitative polarization and orientational analysis 
in different experimental configurations can provide de- 
tailed information on structure, conformation and dy- 
namics of molecules at liquid interfaces [ol. E^ . I23L E3 . l25| . 
The polarization dependent SFG intensity is j22| : 



with 
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in which ni{uji) is the refractive index of phase i at fre- 
quency LUi] /3 is the angle of the laser beam; / is the in- 

(2) 

tensity of the laser beam or the SFG signal. Xeff 

(2) 

effective second order susceptibility; xljk macro- 
scopic susceptibility tensor elements, determined by the 
macroscopic and molecular symmetry |24j : Tj{oji) the lo- 
cal field factor tensor, and e(wi) the unit vector of the 
corresponding optical field, which is responsible for the 
experimental configuration dependence. 

As Wei et al. had demonstrated 9] , when orientational 
motion average is considered, if the orientational motion 
is faster than the vibrational relaxation time scale 1/Tq 
of the qth mode, slow motion average of the orientation 
motion is no longer valid. The slow motion limit gives. 
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and the fast motion gives 



UJ2 — OJq + 



— {D,xD,^Dku) (3) 



(2) 



EE 



a. 



UJ2 



{D,xD,^){Du.) (4) 



in which is the surface density of molecules; aq^x^i, 
and LOq and Tq are the amplitude, resonant frequency and 
damping constant of the qth molecular vibrational mode, 
respectively; and Di^{t) = I ■ ^{t) is the time-dependent 
direction cosine matrix with I — i,j,k for laboratory co- 
ordinates and = X, fijV for molecular coordinates. 

Wei et al. used the step orientational distribution func- 
tion in Eq. as well as other distribution functions, 
such as Gaussian, centered at the surface normal, and 
concluded for a fast motion of the interfacial free OH 
bond in a range of 9m = 51° from the data obtained 
with Vis= 45°, IR=57° p. Otherwise, the ppp and sps 
intensities should be comparable, given the ssp intensity 
is about 10 times of that for ppp in their SFG data. 
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It is clear that such treatment can not explain the much 
stronger sps intensity around ~ 3700cm~^ in Config.4, 
even though it can seemingly explain data in Config.1,2, 
and 3. Using the same parameters and calculation pro- 
cedures as Wei et al, we simulated the ssp, ppp, and sps 
intensities with both slow motion and fast motion limit. 
Fig. [3 shows the results for Config.2 and Config.4. Our 
calculation of Config.2 gives the same results as that by 
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FIG. 2: SFG intensity of the free OH bond simulated with 
both slow motion limit (solid curves) and fast motion limit 
(dotted curves) following the procedure and parameters as 
Wei et al. . In order to compare values in different configu- 
ration, the factor of sec^ /3 is included in the calculation, and 
all intensities are normalized to the ssp intensity in Config.4 
when 9m = 0°. The vertical line indicates the distribution 
width suggested by Wei et al. 



Wei et al. We noticed in the calculation that the 

fast orientational motion centered at interfacial normal 
shall make all SFG intensity in sps spectra vanish for all 
experimental configurations, independent of the range of 
the motion. This is clearly in contradiction to the data 
from Config.4. As Wei et al. concluded, the slow motion 
limit certainly can not explain the relative intensity of the 
3700cTO~^ peak for the ssp, ppp polarizations. Therefore, 
alternative description of the motion and orientation has 
to be invoked. 

Now we consider the case when f{d) is a Gaussian 
function around Oq instead of the interfacial normal. 



fid) 
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Fig. 131 shows simulation of Config.2 and Config.4 as- 
suming a rotationally isotropic interface with cr = 0, us- 
ing the same set of parameters as Wei et al. , and the 
average procedures For data in all four 

configurations in Fig. ^ a consistent 9o = 33° can be 
reached from polarization and orientation analysis |2Cll| . 
This result is in good agreement with the < 38° tilt angle 
reported by Q. Du et al. 0. 

It is easy to show that when ct = 0, fast motion and 
slow motion treatment would converge. Further simu- 
lations show that in order to satisfy the intensities of 



around 3700cto~^ in all four sets of data, a can not ex- 
ceed 15° and Oq changes from 33° to 30° as a increases. 
Within such a relatively small range of cr, simulation re- 
sults for fast and slow motions become indistinguishable. 
Using step function of f{9) centered at 9o does not change 
our general conclusion. This clearly indicates that the as- 
sumption of f(0) centered at interfacial normal is out of 
the question. 

Libration dynamics of the hydrogen bond in liquid wa- 
ter can be as fast as O.lps [l^|22j. Therefore, from the 
observed SFG data our general conclusion is that such 
extremely fast libration dynamics of the tilted free OH 
bond at the air/ water interface is only possible for orien- 
tational motion within a relatively narrow range. 



o 

Li_ 

CO 



^ 3.0- 
w 
c: 
Q) 

.E 
O 

0.0- 



ssp ■ 

ppp\ 

sps ^ 


Config.2 
— 






1 1 1 hn 

PPP\ 

ssp ^ 


— 1 1 1 1 1 1 

Config.4 


sps 









10 20 30 40 50 60 70 80 90 

orientation angtle e (degree) 

FIG. 3; Simulated SFG intensity of the free OH bonds at 
different tilt angle 9o with ct = 0°. All curves are normalized 
to the ssp intensity in Config.4 with 9 — 0°. The vertical 
line indicates the orientation which quantitatively explains 
the observed SFG data. Please note the difference of the 
scales for the two configurations. 



In the SFG literatures on air/water interface, data were 
usually presented for the ssp polarization, except for very 
few cases Polarization analysis were only per- 

formed by Wei et al. Our simulation shows strong de- 
pendence for the ppp intensity of the free OH bond spec- 
tra on the visible incident angle, especially when around 
Vis= 45°. In this range, one or two degree change of 
the visible incident angle can cause significant percent- 
age change of the ppp intensity, as can be indicated from 
data in Config 2 and 3. Wei et a/.'s SFG experiment was 
performed with Vis= 45°. It happened to be that their 
analysis with ultrafast motion in a broad range of orien- 
tation centered around interface normal could just well 
explain this set of data. Therefore, they did not consider 
other possible physical pictures [23. Besides, Richmond 
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et al. suggested two possible qualitative explanations for 
the lack of sps intensity of the free OH stretching mode 
in the SFG measurement [2^ . 

Polarization analysis can further determine the sym- 
metry property of SFG spectral features 0|. Analysis 
indicates that both the sharp peak at ~ 3700cm~^ and 
the broad peak at ~ 3550cm~^ belong to Coov symmetry; 
while both the hydrogen bonded broad peak ~ 3250cm^^ 
and ^ 3450cTO~^ belong to symmetric type of the C2V 
symmetry. Analysis also indicates that the existence of 
interfacial water molecules with two free OH bonds is 
not supported in SFG data Such interfacial water 

molecule shall have C2V symmetry, which should have 
shown clear ssp spectral intensity around 3650cto^^ or 
much stronger ppp spectral intensity at the asymmetric 
stretching mode position of 3750cm~-'^ in at least one of 
the four sets of SFG data 

Thus, ~ 3550cTO~^ peak with C'ocv symmetry, which is 
apparent in the ppp polarization of Config.2 and Config.3, 
can be assigned to the hydrogen bonded OH of the inter- 
facial water molecule with a free OH bond. This is also 
supported with the fact that 3550cm^^ is also the fre- 
quency of the hydrogen bonded OH in water dimer [s^ • 
Orientational analysis of its peak intensities in different 
experimental configurations indicates that this OH bond 
oriented with an angle around ~ 140° from the interface 
normal. Indeed, with the free OH pointing to the vapor 
phase around 35°, the hydrogen bonded OH of the same 
molecule pointing into the liquid phase has to assume an 
orientation around ^ 140°. This is consistent with the 
conclusion that the dipole vector of the interfacial water 
molecule is close to parallel to the interface 0, 0, 0, |^ ■ 
The detail of this analysis shall be presented elsewhere. 

In summary, quantitative analysis of the SFG spectra 
in different polarization and experimental configurations 
for the air/water interface shows that orientational mo- 
tion of the interfacial water molecule is libratory, as fast 
as 0.1 ps it mig ht be mill, only within a limited angu- 
lar range of less than 15° with the tilt angle around 30°. 
Therefore, the air/water interface is quite well-ordered. 
This picture is significantly different from the previous 
conclusion that the interfacial water molecule orientation 
varies over a broad range within the vibrational relax- 
ation time, by Wei et al. using single SFG experimental 
configuration ^ . This progress provided a direct and de- 
tailed physical picture of the orientation and motion for 
the air/water interface. 
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